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Abstract:  Adsorption  capacities  of  adsorbents  are  necessary  for  selecting  and  designing  adsorption  systems 
for  separation  and  removal  processes,  such  as  air  quality  control  devices,  because  they  are  indicators  of  service 
life.  This  paper  describes  the  use  of  the  Dubinin-Radushkevich  (DR)  equation  and  the  quantitative  structure- 
activated  relationship  to  predict  the  equilibrium  adsorption  of  select  organic  vapor  by  activated  carbon  fiber 
(ACF)  adsorbents.  The  DR  isotherm  parameter,  k ,  depends  on  the  adsorbate  as  well  as  the  adsorbent,  and  the 
prediction  for  k  can  be  obtained  indirectly  from  the  affinity  coefficient.  A  correlation  is  developed  to  compute 
the  affinity  coefficient  from  the  modified,  first-order,  valence  molecular  connectivity  index.  This  method  provides 
an  improved  way  to  predict  equilibrium  adsorption  capacities  for  select  volatile  organic  compound  adsorbates 
and  activated  carbon  fiber  adsorbents. 


INTRODUCTION 

Activated  carbon  fiber  (ACF)  in  fabric  or  molded  forms  is 
an  efficient  microporous  adsorbent  for  purifying  indoor  air  and 
separating  organic  vapors  in  many  industrial  gas  streams. 
Compared  with  conventional  granular  or  powdered  activated 
carbon,  ACF  has  higher  adsorption  capacity,  faster  adsorption/ 
desorption  rates,  and  is  suitable  for  in  situ  electrothermal  re¬ 
generation  (Petkovska,  et  al.  1991;  Suzuki  1994;  Rood  et  al. 
1998;  Sullivan  et  al.  1999).  These  advantages  can  lead  to  con¬ 
siderable  reduction  in  the  size  of  adsorption  systems  and  the 
effort  required  to  regenerate  saturated  adsorbent. 

Selection  and  design  of  ACF  adsorption  systems  requires 
knowledge  of  the  adsorption  capacities  for  the  target  organic 
vapors  because  they  are  often  indicators  of  service  life.  Lab¬ 
oratory  measurement  of  adsorption  capacities  may  be  hindered 
by  chemical  toxicity,  equipment  availability,  analytical  cost, 
and  time  requirements  to  complete  the  experiments.  Therefore, 
it  is  advantageous  to  predict  adsorption  capacities  using  an 
adsorption  isotherm  model  with  parameters  based  on  the  phys¬ 
ical  properties  of  the  adsorbent  and  the  adsorbates. 

The  adsorption  capacity  for  an  organic  vapor  is  related  to 
its  partial  pressure  at  equilibrium  by  an  adsorption  isotherm. 
The  Dubinin-Radushkevich  (DR)  equation  is  a  commonly 
used  pure  component  isotherm  model  that  was  developed  from 
the  Polanyi  concept  of  adsorption  potential  and  characteristic 
curves  (Dubinin  et  al.  1991;  Wood  1992;  Dubinin  1985).  The 
DR  equation  is  written  as 
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where  q  =  equilibrium  adsorption  capacity  for  the  adsorbate; 
W0  =  limiting  pore  volume;  p  =  density  of  condensed  liquid 
in  the  micropores;  k  =  DR  isotherm  parameter  for  the  target 
adsorbate;  R  =  universal  gas  constant;  T  =  absolute  tempera¬ 
ture;  P  =  partial  pressure  of  the  adsorbate;  P,  =  adsorbate’s 
bulk  saturation  vapor  pressure;  k,  =  DR  isotherm  parameter 
for  the  reference  adsorbate;  and  (3  =  affinity  coefficient. 

The  isotherm  parameters,  W0  and  k,  must  be  provided  with 
(1)  to  calculate  the  adsorption  capacity  for  an  organic  vapor. 
The  limiting  pore  volume.  W0,  is  known  as  an  adsorbent  struc¬ 
ture  constant.  The  value  of  W0  can  be  measured  using  one 
adsorbate,  and  that  value  can  typically  be  applied  to  the  same 
adsorbent  that  is  adsorbing  other  similar  adsorbates  (Foster  et 
al.  1992;  Qi  et  al.  1998).  Parameter  k  depends  on  the  properties 
of  the  adsorbent  and  adsorbate.  Experimental  determination  of 
k,  or  |3  if  k,  is  known,  can  be  expensive  and  time  consuming, 
and  it  is  impractical  to  measure  these  values  for  all  relevant 
adsorbate-adsorbent  systems.  These  issues  have  caused  interest 
in  developing  methods  to  predict  k  or  p. 

Conventionally,  k  is  assumed  to  be  the  ratio  of  k,  and  [32 
given  by  (2).  Parameter  ks  may  be  determined  using  a  refer¬ 
ence  adsorbate  such  as  benzene,  for  which  (3  =  1.0  (Dubinin 
1975).  The  affinity  coefficient  has  been  correlated  to  the  mo¬ 
lecular  parachor,  the  molar  liquid  volume,  and  the  molar  po¬ 
larizability  of  adsorbates  (Dubinin  1975;  Wood  1992).  Most 
recently,  investigators  (Nirmalakhandan  and  Speece  1993; 
Prakash  et  al.  1994)  proposed  the  quantitative  structure-activ¬ 
ity  relationship  (QSAR)  method  to  predict  k  directly  without 
using  a  reference  adsorbate.  Flowever,  there  is  a  lack  of  con¬ 
vincing  evidence  that  demonstrates  which  of  the  physical  pa¬ 
rameters  mentioned  above  correlate  best  with  (3  and/or  k. 

The  objective  of  this  paper  is  to  investigate  two  QSAR  tech¬ 
niques  for  estimating  the  DR  isotherm  parameter  k  for  adsorp¬ 
tion  of  organic  vapors  by  ACF  adsorbents.  A  correlation  is 
developed  relating  (3  to  the  molecular  connectivity  index.  With 
this  correlation,  k  may  be  calculated  indirectly  using  (2).  Re¬ 
sults  from  the  indirect  QSAR  method  and  the  existing  direct 
QSAR  method  (Nirmalakhandan  and  Speece  1993;  Prakash  et 
al.  1994)  are  then  compared  using  experimental  data  for  ad¬ 
sorption  of  benzene,  acetone,  and  methyl  ethyl  ketone  (MEK) 
onto  three  different  microporous  ACF  adsorbents. 

QSAR  METHODS 

The  modified,  first-order,  valence  molecular  connectivity  in¬ 
dex,  Y,  is  a  molecular  structure  descriptor  originally  devel¬ 
oped  by  Kier  and  Hall  (1986).  As  shown  by  Prakash  et  al. 
(1994),  Y  can  be  calculated  directly  from  the  molecular  struc¬ 
ture  of  a  chemical  by  an  algorithm  requiring  no  experimental 
inputs.  For  this  reason,  Y  was  correlated  to  parameters  k 
and  (3. 
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In  the  direct  QSAR  method,  k  is  assumed  to  be  dependent 
on  the  adsorbate  alone  and  practically  independent  of  the  ad¬ 
sorbent.  By  linear  regression  of  the  best-fit  k  values  versus  Y 
values,  Nirmalakhandan  and  Speece  (1993)  developed  the  fol¬ 
lowing  correlation: 

log  k  =  1.585  -  0.442‘x1',  ( r 2  =  0.924)  (3) 

The  units  of  k  in  (3)  are  10-8  molVcal  '. 

The  best-fit  k  value  was  determined  by  minimizing  the  dif¬ 
ference  between  equilibrium  adsorption  capacities  predicted  by 
the  DR  equation  and  the  experimental  adsorption  isotherm  val¬ 
ues  for  select  combinations  of  organic  vapors  and  granular  and 
powdered  activated  carbon  adsorbents.  As  pointed  out  by 
Prakash  et  al.  (1994),  the  advantage  of  this  approach  is  that 
no  reference  chemical  and  experimental  inputs  are  required  to 
estimate  k.  However,  to  use  the  DR  Equation,  W0  must  be 
known.  The  value  for  W0  may  be  determined  from  isotherm 
data  using  a  representative  adsorbate  such  as  benzene  or  water 
vapors  (Foster  et  al.  1992;  Qi  et  al.  1998). 

The  indirect  QSAR  method  uses  the  conventional  approach 
that  considers  interactions  between  the  adsorbate  and  adsor¬ 
bent  described  by  (2)  with  the  independent  parameters,  k,  and 
p.  The  following  correlation  relates  (3  to  Y: 

P  =  0.3690Y  +  0.2729,  ( r 2  =  0.967)  (4) 

Eq.  (4)  was  developed  by  linear  regression  of  the  best-fit  k 
values  versus  Y  values,  where  the  best-fit  k  values  were  de¬ 
termined  using  the  experimental  data  of  Cal  (1995)  for  ad¬ 
sorption  of  benzene,  acetone,  and  MEK  onto  three  ACF  ad¬ 
sorbents.  The  adsorbents  were  manufactured  from  cross-linked 
phenol-aldehyde  fibers  by  American  Kynol,  Inc.  (Table  1).  The 
limiting  pore  volume,  W0,  was  determined  simultaneously  with 
ks  using  benzene  as  a  reference  chemical. 

ANALYSIS  OF  QSAR  METHODS 

The  two  QSAR  methods  were  evaluated  using  the  experi¬ 
mental  data  of  Cal  (1995)  for  adsorption  of  benzene,  acetone, 
and  MEK  onto  three  different  ACF  adsorbents  (Table  1).  The 
adsorption  experiments  were  performed  using  a  gravimetric 


balance  (Cahn  Model  2000)  at  room  temperature  (~25°C)  and 
a  total  pressure  of  1  atm  (Cal  1995).  There  are  nine  adsorption 
isotherms  that  were  generated  with  more  than  100  data  points 
using  organic  vapor  concentrations  in  the  ppmv  range.  The 
saturated  vapor  pressures  of  the  organic  compounds  were  es¬ 
timated  using  the  Wagner  equation  (Reid  et  al.  1986),  and  the 
density  of  the  organic  compound  as  a  bulk  liquid  was  deter¬ 
mined  by  Beaton  and  Hewit  (1989).  For  each  adsorbent,  W0 
and  ks  (with  p  =  1.0)  were  determined  by  minimizing  the  dif¬ 
ference  between  the  results  from  the  DR  equation  and  the  ben¬ 
zene  adsorption  isotherm  data. 

Values  for  k  estimated  with  the  direct  QSAR  method  (3)  are 
compared  with  the  corresponding  k  values  estimated  with  the 
best-fit  k  values  in  Table  2.  For  adsorption  of  benzene,  acetone, 
and  MEK  onto  ACF-20  and  ACF-25,  the  estimates  by  the  di¬ 
rect  QSAR  and  the  best-fit  k  methods  agree  within  a  maximum 
difference  of  16%.  However,  the  differences  in  estimated  k 
values  are  much  higher  for  the  three  organic  vapors  that  ad¬ 
sorbed  onto  ACF-15  (63-96%  difference).  This  may  result 
from  not  considering  the  adsorbate-adsorbent  interactions. 
Consequently,  the  DR  equation  with  k  calculated  by  the  direct 
QSAR  method  significantly  underestimates  the  adsorption  iso¬ 
therms  for  the  organic  vapor-ACF-15  systems  (Fig.  1).  The 
average  errors  are  27,  35,  and  62%  for  benzene,  MEK,  and 
acetone,  respectively,  when  using  the  direct  QSAR  method  to 
estimate  k. 

Table  2  also  compares  k  values  estimated  with  the  indirect 
QSAR  method  and  the  best-fit  k  values.  The  k  values  estimated 
indirectly  using  (2)  and  (4)  appear  to  agree  well  with  the  best- 
fit  k  values  for  the  three  adsorbates  onto  the  three  ACF  ad¬ 
sorbents  (11.  5,  and  2%  maximum  differences  for  ACF-15, 
ACF-20,  and  ACF-25,  respectively).  This  may  be  attributed  to 
the  consideration  of  both  adsorbent  and  adsorbate  by  the  in¬ 
direct  QSAR  method.  With  the  estimated  k  values  and  the  DR 
equation,  the  adsorption  isotherms  are  computed  and  com¬ 
pared  with  the  measured  isotherm  data  in  Figs.  2-4.  The  DR 
equation  agrees  well  with  the  measured  isotherm  data  for  ad¬ 
sorption  of  benzene  and  acetone  onto  the  three  ACF  adsor¬ 
bents.  Similar  results  are  shown  for  MEK  on  ACF-15  and 
ACF-20.  However,  the  DR  equation  underestimates  the  ad- 


TABLE  1 .  Properties  of  Activated  Carbon  Fiber  Adsorbents 


Property 

(1) 

ACF-15 

(2) 

ACF-20 

(3) 

ACF-25 

(4) 

BET  specific  surface  area  (P/P,  for  N,)  =  10  —  5  to  0.4  at  77  K)  (m2/g)“ 

730 

1,330 

1.860 

C/H/O/N  (mass  percent)11 

92.8/1.04/6.12/<0.05 

95.4/.68/3.92/<0.05 

95. 4/0.59/3. 07/<0.05 

Pore  size  with  maximum  volume  (n  m )' 

7.0 

— 

9.0 

“Cal  1995. 
bFoster  et  al.  1992. 
cSun  et  al.  1998. 


TABLE  2.  Comparison  of  Computed  and  Best-Fit  DR  Isotherm  k  Values3 


Adsorbent 

(1) 

Chemical 

(2) 

w0 

(cm3/g) 

(3) 

V 

(4) 

3 

(5) 

Best-fit 

k 

(mole/KJ)2 
(k  x  1 03) 

(6) 

Direct  QSAR 
k 

(mole/KJ)2 
(k  x  1 03) 

(7) 

Indirect  QSAR 
k 

(mole/KJ)2 

(8) 

ACF-15 

Benzene 

0.457 

2.0 

1.01 

1.77“ 

2.87 

1.73  X  10“3 

Acetone 

0.457 

1.2 

0.72 

3.85 

6.47 

3.46  X  10“3 

MEK 

0.457 

1.76 

0.92 

1.87 

3.66 

2.08  X  10“3 

ACF-20 

Benzene 

0.715 

2.0 

1.01 

3.20 

2.87 

3.13  X  10“3 

Acetone 

0.715 

1.2 

0.72 

5.96 

6.47 

6.24  X  10“3 

MEK 

0.715 

1.76 

0.92 

3.70 

3.66 

3.74  X  10“3 

ACF-25 

Benzene 

0.563 

2.0 

1.01 

3.42 

2.87 

3.35  X  10“3 

Acetone 

0.563 

1.2 

0.72 

6.57 

6.47 

6.67  X  10“3 

MEK 

0.563 

1.76 

0.92 

4.08 

3.66 

4.02  X  10“3 

* k  =  1.77  X  10“3  (mole/KJ)2. 
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sorption  capacity  for  MEK  on  ACF-25  as  the  vapor  concen¬ 
tration  increases.  The  largest  average  error  for  any  of  the  com¬ 
puted  isotherms  in  Figs.  2-4  is  14%. 

It  is  important  to  note  that  the  original  direct  QSAR  method 
was  developed  using  a  broad  database  while  the  indirect 
QSAR  method  presented  here  was  developed  using  the  limited 
ACF  data  of  Cal  (1995)  to  which  it  has  also  been  compared. 
This  may  have  contributed  to  the  larger  regression  coefficients, 
r2,  of  the  indirect  method  compared  with  the  direct  method. 

In  this  study,  benzene  was  used  as  the  reference  chemical, 
although  theoretically  any  chemical  can  be  chosen.  Benzene 
is  nonpolar,  while  acetone  and  MEK  are  polar.  The  reasonable 
approximations  of  the  adsorption  data  for  these  three  organic 
vapors  by  the  indirect  QSAR  method  indicate  that  the  chem¬ 
ical  polarity  has  a  limited  effect  on  p. 

Table  3  compares  the  values  for  (3  estimated  using  three 
different  predictive  methods  for  selected  organic  vapors.  These 
methods  are  correlations  relating  [3  to  the  molecular  connec¬ 
tivity  index,  the  molecular  parachor,  and  the  molar  polariza¬ 
bility,  respectively.  A  detailed  description  of  the  molecular 
parachor  method  and  the  polarizability  method  are  available 
in  the  literature  (Dubinin  1975;  Dubinin  1991;  Cal  1995).  As 
shown  in  Table  3,  the  three  methods  provide  similar  estimates 
for  (3.  The  estimated  [3  values  by  the  indirect  QSAR  method 
can  be  linearly  correlated  with  those  estimated  by  the  other 
two  methods  with  a  regression  coefficient  of  r2  =  0.97.  This 


FIG.  1.  Measured  and  Computed  Adsorption  Isotherms  Using 
DR  Equation  with  k  Values  from  Direct  QSAR  Method  for  Ben¬ 
zene,  Acetone,  and  MEK  on  ACF-15  Systems 


0  200  400  600  800  1000 


Vapor  Concentration,  ppmv 

FIG.  2.  Measured  and  Computed  Adsorption  Isotherms  Using 
DR  Equation  with  /(Values  from  Indirect  QSAR  Method  for  Ben¬ 
zene  on  3  ACF  Adsorbents 
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FIG.  3.  Measured  and  Computed  Adsorption  Isotherms  Using 
DR  Equation  with  k  Values  from  Indirect  QSAR  Method  for  Ace¬ 
tone  on  3  ACF  Adsorbents 


FIG.  4.  Measured  and  Computed  Adsorption  Isotherms  Using 
DR  Equation  with  /(Values  from  Indirect  QSAR  Method  for  MEK 
on  3  ACF  Adsorbents 


TABLE  3.  Values3  for  Affinity  Coefficient  Estimated  Using  Se¬ 
lect  Methods  and  Organic  Vapors 


Chemical 

(1) 

Indirect  QSAR 

P 

(2) 

Parachor 

P 

(3) 

Polarizability 

P 

(4) 

Acetaldehyde 

0.53 

0.65 

0.43 

Acetone 

0.72 

0.78 

0.61 

Benzene 

1.01 

1.00 

1.00 

Ethyl  Chloride 

0.83 

0.72 

0.873 

MEK 

0.92 

0.96 

0.81 

Tetrachloride 

1.11 

— 

1.12 

Toluene 

1.16 

1.19 

1.20 

P-xylene 

1.31 

1.38 

1.38 

“Based  on  benzene  as  reference  chemical. 


suggests  that  the  indirect  QSAR  method  should  be  useful  be¬ 
yond  the  limited  database  tested  here. 

SUMMARY  AND  CONCLUSIONS 

This  paper  demonstrates  the  use  of  the  Dubinin-Radushke- 
vich  (DR)  equation  and  the  quantitative  structure-activated  re¬ 
lationship  (QSAR)  methods  to  estimate  the  adsorption  capac¬ 
ities  of  select  organic  vapors  and  activated  carbon  fiber  (ACF) 
adsorbents.  Two  QSAR  methods  are  evaluated  using  the  ex¬ 
perimental  data  for  adsorption  of  benzene,  acetone,  and  methyl 
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ethyl  ketone  onto  three  different  ACF  adsorbents.  It  is  sug¬ 
gested  that  the  DR  isotherm  parameter  k  depends  on  the  ad¬ 
sorbate  as  well  as  the  adsorbent,  and  an  improved  prediction 
for  k  may  be  obtained  indirectly  from  the  affinity  coefficient. 
A  correlation  is  presented  to  compute  the  affinity  coefficient 
from  the  modified,  first-order,  valence  molecular  connectivity 
index. 
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APPENDIX  II.  NOTATION 

The  following  symbols  are  used  in  this  paper: 

k  =  DR  isotherm  parameter  for  target  adsorbate  (mole/kj)2; 
k,  =  DR  isotherm  parameter  for  reference  adsorbate; 

P  =  partial  pressure  of  target  adsorbate; 

P0  =  saturation  partial  pressure  of  target  adsorbate; 
q  =  adsorption  capacity  (mg/g); 

R  =  universal  gas  constant; 
r2  =  regression  coefficient; 

T  =  absolute  temperature  (K); 

W0  =  limiting  pore  volume  (cmVg); 

(3  =  affinity  coefficient; 

p  =  density  of  condensed  liquid  in  micropores  (g/cnv);  and 
Y  =  molecular  connectivity  index. 
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